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ABSTRACT 

Aims. PSR J 1846-0258 is a young rotation-powered pulsar with one of the highest surface magnetic field strengths, located in the 
centre of SN-remnant Kes-75. In June 2006 a magnetar-like outburst took place. Using multi-year RXTE and INTEGRAL observa- 
tions covering also the epoch of the outburst we aim to study the temporal and spectral characterisitics of PSR J1846-0258 over a 
broad ~ 3-300 keV energy range to derive constraints on theoretical scenarios aiming to explain this schizophrenic behaviour 
Methods. We explored all publically available RXTE observations of PSR J1846-0258 to generate accurate ephemerides over the 
period January 30, 2000 - November 7, 2007. Phase folding procedures yielded pulse profiles for RXTE PCA (~ 3 - 30 keV), RXTE 
HEXTE (~ 15 - 250 keV) and INTEGRAL ISGRI (~ 20 - 300 keV). The pulsed spectrum over the full - 3-300 keV energy range 
was derived, as well as the total spectrum (including the Pulsar Wind Nebula) over the 20-300 keV band with the ISGRI. The timing, 
spatial and spectral analyses have been applied for epochs before, during and after the magnetar-like outburst to study the evolution 
of the high-energy characteristics. 

Results. ISGRI detected PSR J1846-0258/Kes-75 before outburst during 2003-2006 with a power-law-shape spectrum over the 20- 
300 keV energy range with photon index T = 1.80 ± 0.06 and energy flux (20-300 keV) of (6.62 ± 0.35) x 10"" erg/em's. More 
than 90 days after the onset of the outburst, still during the decay phase, the same spectral shape was measured (T = l.TS'^g'll) with 
an indication for a 52% (2.3(t) enhanced total emission, while one year after the outburst the hard X-ray non-thermal emission of 
PSRJ1846-0258/Kes-75 was found to be back to its pre-outburst values. PCA monitoring of PSR J1846-0258 before the outburst 
yielded phase-coherent ephemerides confirming the earlier derived breaking index of the spindown. During the outburst incoherent 
solutions have been derived. We showed that the radiative outburst was triggered by a major spin-up glitch near MJD 53883 ± 3 with 
a glitch size Av/v in the range (2.0 - 4.4) x 10"*. Using all pre-outburst observations of ISGRI and HEXTE for the first time pulse 
profiles have been obtained up to 150 keV with a broad single asymmetric pulse. The pulse shape did not vary with energy over the 
2.9-150 keV energy range, nor did it change during the magnetar-like outburst. The time-averaged pre-outburst ~ 3-300 keV pulsed 
spectrum measured with the PCA, HEXTE and ISGRI was fitted with a power-law model with F = 1.20 ± 0.01. A fit with a curved 
power-law model gives an improved fit. Around 150 keV the pulsed fraction approaches 100%. The first 32 days of the magnetar-like 
outburst the 3-30 keV pulsed spectrum can be represented with two power laws, a soft component with index T, = 2.96 ± 0.06 and a 
hard component with the pre-outburst value T;, ~ 1.2. Above ~ 9 keV all spectra during outburst are consistent with the latter single 
power-law shape with index ~ 1.2. The 2-10 keV flux increased by a factor ~ 5 and the 10-30 keV flux increased with only 35%. 
After ~ 120 days the soft outburst and the enhancement of the hard non-thermal component both vanish. 

Conclusions. The varying temporal and spectral characteristics of PSR J1846-0258 can be explained in a scenario of a young high-B- 
field pulsar in which a major glitch triggered a sudden release of energy. Resonant cyclotron upscattering could subsequently generate 
the decaying / cooling soft pulsed component measured during outburst between 3 and 10 keV. The (variation in the) non- thermal 
hard X-ray component can be explained with synchrotron emission in a slot-gap or outer-gap pulsar model. 

Key words. Stars: neutron - pulsars: individual PSR J1846-0258, PSR B 1509-58, 4U 0142+61, IRXS J 170849-4009 10 - X-rays: 
general - Gamma rays: observations 



1. Introduction 

PSR J 1846-0258 was discovered in a timin g analysis of X- 
ray data from RXTE and ASCA by Gotthel f et all (l2000l) . It 
is a relatively slow rotation-powered pulsar, P ~ 324 ms, and 
has one of the highest surface magnetic field strengths, B ~ 
4.9 X lO'-' G, assuming standard magnetic dipole breaking, just 
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above the quantum critical field strength B^r = nigr' jefi of 
4.413 X 10'"* G. Furthermore, it has the smallest characteristic 
age, T ~ 723 y, of all known pulsars and a spin-down lumi- 
nosity (Esd) 8.2 X 10"^^ er g/s. At radio frequencies it remains 
undetected ([Archibald et all 2008). PSR J1846-0258 is located 
at t he centre of Supernova Remnant (SNR) Kes 75 (G29.7-0.3; 
,Kesteveni[T968h and was resolved in a Chandra observation as 
a bright X-ray source surrounded by a diffuse pulsar wind neb- 
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Table 1. INTEGRAL observation characteristics for PSR J 1846-0258 



Revs. 


Date begin 


Date End 


MJD GXr exposure 


Eff.'' exposure 


# Sew-' 










(Ms) 


(Ms) 










Pre-outburst 








049-070 


10-03-2003 


13-05-2003 


52708-52772 


1.3737 


0.6954 


657 


109-123 


04-09-2003 


18-10-2003 


52886-52930 


0.3434 


0.1869 


180 


172-233 


11-03-2004 


10-09-2004 


53075-53258 


1.1416 


0.5047 


402 


236-250 


18-09-2004 


01-11-2004 


53266-53310 


0.6488 


0.3142 


343 


251-315 


02-11-2004 


14-05-2005 


53311-53504 


0.6681 


0.3945 


315 


345-379 


10-08-2005 


21-11-2005 


53592-53695 


1.0405 


0.4418 


347 


40/-44i 


12-02-2006 


26-05-2006 


53778-53881 


0.7707 


0.4604 


355 


049-441 


10-03-2003 


26-05-2006 


52708-53881 


5.9868 


2.9979 


2599 








Outburst/Post-outburst 






AHA "^ni 


Jl-US-zUUo 


1 Q 11 ic\r\(^ 
IV- i l-ZUUD 


53978-54058 


0.3721 




1 /I G 


537-561 


07-03-2007 


20-05-2007 


54166-54241 


0.3908 


0.1561 


130 


592-603 


19-08-2007 


23-09-2007 


54331-54366 


0.6011 


0.2587 


168 


474-603 


31-08-2006 


23-09-2007 


53978-54366 


1.3640 


0.6283 


446 








All observations 








049-603 


10-03-2003 


23-09-2007 


52708-54366 


7.3508 


3.6262 


3045 



' Total Good-Time-Interval exposure of the used observations 

^ Effective exposure on PSR J 1846-0258 corrected for off-axis sensitivity reduction 

' Number of used Science Windows, see text 



ula (PWN; iHelfand et all |2003|) . The distance of SNR Kes 75 
and PSR J1846-0258 has been under discu ssion for a long time. 
The most recent estimates are bv lLeahv & T ian (2008), 5. 1 to 7.5 
kpc, and lSu et"an (12008 ^, ~ 10.6 kpc. For a discussion on earlier 
estimates we refer to their papers. In this work we adopt a dis- 
tance of 10 kpc, and write the distance to Kes 15 as d - dio x 10 
kpc. 

Based on analysis of more than five years of RXTE moni- 
toring data it appeared that PSR J 1846-0258 behaves as a very 
sta ble rotator with a breaking index of « = 2.65 + 0.01 (v oc 
v"; [Livingstone et all [2006). INTEGRAL detected hard X-ray 
point-source emission from the PWN plus PSR J1846-0258 up 
to ~ 200 keV (McB ride et al.. 2008). and HESS discovered the 
PWN in the TeV band loiannati-Atai et alll2007l) . 

Recently, Ku mar & Safi-Harbl ( l2008l) reported the detection 
of a dramatic brightening of the pulsar during June, 7-12, 2006 in 
Chandra observations of Kes 75. The pulsar's spectrum (pulsed 
plus DC emission) softened considerably from a power-law in- 
dex r ~ 1.32 measured i n 2000 to 1.97 during the outburst 
in 2006. Gavriil et al. (2008) used RXTE observations to show 
that this prompt radiative event exhibited a recovery from the 
pulsed-flux enhancement which can be modeled as an exponen- 
tial decay with 1/e time constant of 55.5 + 5.7 day and a to- 
tal 2-60 keV energy rele ase of (3.8 - 4 . 8) x 10"^' erg, adopt- 
ing a distance of 6 kpc dLeahy & TianL 12008 ) and assuming 
isotropic emission. For the revised distance estimate of ~ 10 
kpc the released energy becomes (1.1 - 1.3) x IQ^'d^^^. Such 
an energy release is of similar magnitude as has been detected 
from a few Anomalous X-ray Pulsars (AXPs). The outburst was 
accompanied by an unprecedente d change in timing behaviour 
Furthermore. iGavriil et al.l (l2008h also discovered short (< 0.1 
s) magnetar-like bursts during the outburst, 4 at the beginning 
and 1 near the end. This otherwise very stably behaving rotat- 
ing neutron star, both temporally and spectrally, switched sud- 
denly to unpredictable magnetar-like behaviour This is the first 
time that such dual characteristics have been seen. Therefore, the 



magnetar-like outburst of PSR J1846-0258 warrants a more de- 
tailed study to obtain tighter constraints for modelling this tem- 
porary magnetar-like manifestation of a rotation-powered pulsar. 



The main motivation of this paper is to study the X-ray 
spectral and timing characteristics before, during and after the 
outburst including for the first time the hard X-ray regime 10- 
300 keV. In Section 2 we introduce the instruments used in this 
work, the PCA and HEXTE aboard NASA's Rossi X-ray Timing 
Explorer (RXTE) and IBIS-ISGRI aboard ESA's INTErnational 
Gamma-Ray Astrophysics Laboratory (INTEGRAL). In Section 
3 we present an INTEGRAL deep-exposure sky map of the 
Scutum region (20-70 keV) reveahng PSR J 1846-0258 /Kes 75 
surrounded by several near-by sources, and derive for ener- 
gies above 20 keV three spectra for the total emission, PWN 
and PSR J 1846-0258, namely, pre-outburst, during outburst and 
post-outburst. Section 4 deals with the timing studies. First, the 
high statistics in the PCA data are exploited to derive phase- 
coherent ephemerides for PSR J 1846-0258, pre-outburst, for 
part of the outburst and post-outburst, as well as incoherent so- 
lutions during the outburst, revealing a major spin-up glitch at 
the start of the outburst. These (in)coherent timing solutions 
are subsequently used to derive pulse profiles over the energy 
windows covered by the PCA, HEXTE and ISGRI for the pre- 
outburst observations, and using just the PCA for the outburst 
period. In Section 5 we derive the spectrum of the pre-outburst 
total-pulsed emission using PCA, HEXTE and ISGRI data and 
compare this with the total, PWN and PSR J 1846-0258, spectra 
measured with Chandra and INTEGRAL-ISGRI. Furthermore, 
we use PCA data to study the evolution of the spectrum of the 
pulsed signal pre-outburst, during the outburst and post outburst. 
Finally, the results are summarized in Section 6 and discussed in 
Section 7. 
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2. Instruments 

2.1. RXTE 

In this study extensive use is made of data from monitoring 
observations of PSR J1846-0258 with the two non-imaging X- 
ray instruments aboard RXTE, the Proportional Counter Array 
(PCA; 2-60 keV) and the High Ener gy X-ray Tirning Experiment 
(HEXTE; 15-250 keV). The PCA djahoda et all [1996) consists 
of five collimated xenon proportional counter units (PCUs) with 
a total effective area of ~ 6500 cm^ over a ~ 1° (FWHM) field of 
view. Each PCU has a front Propane anti-coincidence layer and 
three Xenon layers which provide the basic scientific data, and 
is sensitive to photons with energies in the range 2-60 keV. The 
energy resolution is about 18% at 6 keV. All data used in this 
work have been collected from observations in GoodXenon or 
GoodXenonwithPropane mode allowing high time resolution 
(0.9/zs) studies in 256 spectral channels. 

The HEXTE instrument dRothschild et all 1 1998b consists of 
two independent detector clusters, each containing four Na(Tl)/ 
CsI(Na) scintillation detectors. The HEXTE detectors are me- 
chanically collimated to a ~ 1 ° (FWHM) field of view and cover 
the 15-250 keV energy range with an energy resolution of ~ 15% 
at 60 keV. The collecting area is 1400 cm^ taking into account 
the loss of the spectral capabilities of one of the detectors. The 
best time resolution of the tagged events is 7.6//S. In its default 
operation mode the field of view of each cluster is switched on 
and off source to provide instantaneous background measure- 
ments. Due to the co-alignment of HEXTE and the PCA, they 
simultaneously observe the same field of view. 

2.2. INTEGRAL 

The INTEGRAL spacecraft dWinkler et al.L |2003|) . launched 
17 October 2002, carries two main y-ray instruments: a 
high-angular-resolution imager IBIS (Ubertini et al., 2 0031) and 
a hig h-energy-resolution spectrometer SPI jVedrenne et al.L 
l2003h . Both instruments are equiped with coded aperture masks 
enabling image reconstruction in the hard X-ray/soft y-ray band. 

Driven by sensitivity considerations, we used only data from 
the IN TEGRAL Soft Gamma-Ray Imager ISGRI (iLebrun et al.L 
l2003h . the upper detector layer of IBIS, sensitive to photons with 
energies in the range ~20 keV - 1 MeV. With an angular reso- 
lution of about 12' and a source location accuracy of better than 
1' (for a > lOcr source) this instrument is able to locate and sep- 
arate high-energy sources in crowded fields within its 19° x 19° 
field of view (50% partially coded) with an unprecedented sen- 
sitivity (~ 960 cm^ at 50 keV). Its energy resolution of about 
7% at 100 keV is amply sufficient to determine the (continuum) 
spectral properties of hard X-ray sources in the ~ 20 - 300 keV 
energy band. 

The timing accuracy of the ISGRI time stamps recorded on 
board is about 61yL(s. The time alignment b etween INTEGRAL 
and RXTE is better than ~ 25fis (see e.g. iKuiper et all 120031: 
iFalanga et all l2005h . Given the fact that th e accuracy of th e 
RXTE clock in absolute time is about 2/js ("Rots et all , |2004|) . 
this implies that the INTEGRAL absolute timing is better than 
~ 27yus. Data from regular INTEGRAL Crab monitoring obser- 
vations show that the clock behaviour is stable, allowing timing 
studies of weak pulsars by accumulating data taken over many 
years. 

In its default operation mode INTEGRAL observes the sky 
in a dither pattern with 2° steps, which could be rectangular e.g. 
a 5 X 5 dither pattern with 25 grid points, or hexagonal with 



7 grid points (target in the middle). Typical integration times 
for each grid point (pointing/sub-observation) are in the range 
1800 - 3600 seconds. This strategy drives the structure of the 
INTEGRAL data archive which is organised in so-called sci- 
ence windows (Sew) per INTEGRAL orbital revolution (lasting 
for about 3 days) containing the data from all instruments for a 
given pointing. 

3. ISGRI Imaging and Total Spectra 

3.1. Deep ISGRI Sky Map of Scutum Region 

In this work we exploited the good imaging capability of ISGRI, 
which allowed us to separate PSR J1846-0258/Kes 75 from 
nearby sources in this crowded region in the first Galactic quad- 
rant. ISGRI sky mosaics (a combination of deconvolved images 
of single science windows) were produced in 10 logarithmically 
binned energy bands covering the 2 0-300 keV energy wind ow. 
We used the imaging software tools jGoldwurm et all l2003h of 
the Offline Scientific Analysis (OS A) package version 5.1 dis- 
tributed by the IN TEGRAL Science Data Centre (ISDC; see e.g. 
ICourvoisier et all 120031) . The end product of this analysis pro- 
vides the count rate, its variance, exposure and significance for 
all 10 energy bands over the (deconvolved) mosaiced sky field. 

In the observation selection we only accepted science 
windows with instrument pointings within 14°5 from the 
PSR J 1846-0258 position. This ensures that (a part of) the de- 
tector plane is illuminated by the target. The resulting list is 
further screened on erratic count rate variations, indicative for 
particle-induced effects due to Earth-radiation-belt passages or 
solar-flare activities, by inspecting visually the count rate in 
20-30 keV band versus time. Science windows showing erratic 
count-rate variations are excluded for further analysis. In the 
event-selection pro cess we only use ev ents with rise times be- 
tween 7 and 90 (see lLebrun et alll2003L for definition), detected 
in non-noisy ISGRI detector pixels. 

The INTEGRAL observations used in this work are listed 
in Table [1] Combining all these publicly available data from 
INTEGRAL revolutions 49 up to and including 603 we ob- 
tained a total exposure of 7.35 Ms (Good-Time-Intervals, GTI). 
Fig. [1] shows for the 20-70 keV range (weighted sum of 5 
energy bands) the significance image of the field surrounding 
PSR J1846-0258/Kes 75. In this image we adopted a detection 
significance threshold of 4cr. Clearly visible at the center of 
the image with high significance (17.4cr) is the emission from 
PSRJ1846-0258 and its PWN (not resolved). Other interest- 
ing high-energy sources are also indicated in this 10° x 5° im- 
age of the Scutum field, particularly Anomalous X-ray Pulsar 
IE 1841-045, embedded in SNR Kes 73, and AX J1838.0- 
06 55, which harbours a young energetic pulsar as found recently 
by Gotthelf & HalpernI (1200 8). The identifications of the other 
sources are given in the figure caption. 

The total spectrum of PSR J1846-0258 and Kes 75 can be 
constructed from the 10 count-rate and variance maps by ex- 
tracting the (dead-time corrected) rates and uncertainties at the 
location of PSR J 1846-0258. These values are normalized to the 
count rates measured for the total (nebula and pulsar) emission 
from the Crab in similar energy bands using Crab calibration 
observations during INTEGRAL revolutions 102 and 103. From 
the ratios and the photon spectrum of the total emission from the 
Crab, we can derive the total high-energy photon spectrum of 
the PWN plus PSR J 1846-0258 (pulsed and any unpulsed point 
source component) without detailed knowledge of the ISGRI 
energy response. For the total Crab photon spectrum we use 
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Fig. 1. A 20-70 keV image in Galactic coordinates of the Scutum region centered on (l,b)=(30,0) combining IBIS ISGRI data from 
INTEGRAL revolutions 49 up to and including 603 (March 10, 2003 - Sept. 23, 2007). The total GTI exposure amounts 7.35 Ms 
for this mosaic image, while the effective exposure on PSR J1846-0258 equals 3.63 Ms. Only source features with significances 
above 4cr are shown. The white circle indicates PSR J1846-0258/Kes 75 . The yellow encircled sources are AXP IE 1841-04 5 (see 
iKuiperet all |2004 l2006l) and the rotation-powered pulsai- AX J1838.0-0655 jGotthelf & Halpernl l2008t iKuiper et all l2008l) . The 
other sources indicated by numbers are: [1] - GS 1843+009,[2] - IGR J18490-0000, [3] - IGR J18485-0047, [4] - XTE J1855-026, 
[5] - 3A 1845-024, [6] - IGR J 18462-0223, [7] - IGR J 18483-03 11, [8] - IGR J 18450-0435 and [9] - AX J1841.0-0535 



the b roken-power-law spectrum derived by iJourdain & Roque"sl 
(I2OO8) based on INTEGRAL-SPI observations of the Crab at 
energies between 23 and 1000 keV. The latest Crab cross calibra- 
tions between SPI and IBIS-ISGRJ provided consistent results. 

3.2. Total pre-outburst spectrum in 20-300 keV band 

Following the above procedure, the pre-outburst (Revs. 49-441, 
see Table [U IBIS-ISGRI spectrum was derived and can be fit- 
ted with a power-law model over the 20-300 keV energy range 
yielding a photon index F = 1.80 + 0.06 and an energy flux 
(20-300 keV) of (6.62 + 0.35) x 10"" erg/cm^s (xl = 0.878 
for 10-2 degrees of freedom; errors are Icr, throughout this pa- 
per). This index is a slightl y harder, but consiste nt with the 
value 2.0 + 0.2 obtained bv iMcBride et all (l2008h . who used 
less data. The 20-100 keV energy flux derived in this work 
(3.47+0.19)xl0"" erg/cm^s, however, is about 20% hig her flian 
the val ue, (2.9+y j) x 10"'' erg/cm^s, obtained bv McBride et all 
( I2OO8I) . This discrepancy can be fully attributed to their use of 
obsolete energy response matrices for IBIS ISGRI in the OSA 
5.1 environment, which yield inconsistent spectral results for the 
total Crab spectrum. Our ISGRI flux measurements for the total 
pre-outburst emission between 20 and 300 keV from the PWN 
plus PSR J1846-0258 are shown in Fig.|6]as purple data points 
with the best power-law fit superposed, and will be discussed 
later together with the spectrum of the pulsed emission. 

3.3. Total spectrum (20-300 keV) during outburst decay 

Unfortunately, INTEGRAL was not pointing to the Scutum re- 
gion during the early phase of the outburst. The first exposure 



started on MJD 53978 (Rev. 474), more than 90 days after the 
onset of the outburst, with the source still showing enhanced 
pulsed flux levels about 40% higher than the pre-outburst level 
(see e.g. top panel of Fig. |7]i. A spectrum of the total emission 
from PSRJ1846-0258/Kes 75 from INTEGRAL observations 
taken between revolutions 474 and 501 (0.2135 Ms effective on- 
source exposure time; more than 95% of the exposure has been 
accumulated during Revs. 474 & 475) could adequately be de- 
scribed by a power-law with photon index F - 1.75^y 3j and 
energy flux (20-100 keV) of (5.27+JJ-^) x lO"" erg/cm^s. The 
spectral index is fully consistent with that of the pre-outburst 
spectrum, but the total flux is about 52% higher than the pre- 
outburst flux value, in line with the increased pulsed flux level 
shown in PCA data at the same epoch at lower energies. The in- 
crease, however, represents only a ~ 2.3cr eff'ect relative to its 
pre-outburst value. Therefore, we can not claim a significant in- 
crease of the total 20-100 keV flux, but only an indication for 
enhanced total emission in the 20-100 keV band. 



3.4. Total post-outburst spectrum in 20-300 keV band 



About one year after the onset of the outburst INTEGRAL ob- 
served PSR J1846-0258/Kes 75 for an effective exposure of 
0.415 Ms during MJD 54166-54366 (INTEGRAL Revs. 537- 
603, see Table[TJ. Its derived total 20-300 keV spectrum could be 
described by a power-law model with F = 1 -90^0 35 and a 20-100 
keV energy flux of (3.64;^q 5^) x 10"" erg/cm^s, fully consistent 
with the pre-outburst measurement. 
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Entry 
# 


Start 
[MJD] 


End 
[MJD] 


to. Epoch 
[MJD,TDB] 


V 

[Hz] 


V 

xlO-" Hz/s 


V 

xlO^^' Hz/s- 




Validity range 
(days) 










Pre-Outburst 













51286 


51290 


51286.0 


3.08273789(23) 


-6.98(15) 


0.0 (fixed) 


0.6850 


4 


li 


51574 


52237 


51574.0 


3 081061350908) 


-6 73176(1) 


3.821(4) 


0.9825 


663 


2^ 


52369 


52837 


52523.0 


3.0755528807(5) 


-6.707486(9) 


3.86(3) 


0.6501 


468 


3 


52915 


53148 


52915.0 


3.0732836450(59) 


-6 6963n) 


4.98(13) 


0.1971 


233 


4 


53030 


53465 


53030.0 


3.0726185382(34) 


-6.69068(3) 


3.93(2) 


0.5030 


435 


53 


53464 


53880 


53464.0 


3.0701124555(45) 


-6.67599(4) 


3.78(3) 


0.8985 


416 










Outburst-B 










6 


53978 


54033 


53997.0 


3.0669883665(67) 


-6.8688(7) 


161(15) 


0.3092 


55 








Post-Outburst; coherence fully 


caught-up 








7 


54126 


54229 


54126.0 


3.066215110(19) 


-6.8726(8) 


202(2) 


0.2658 


103 


8 


54237 


54340 


54237.0 


3.065565121(15) 


-6.7297(7) 


15(2) 


0.9120 


103 


9 


54340 


54411 


54340.0 


3.064967085(26) 


-6.707(2) 


49(6) 


0.8056 


71 



' Small glitch reported by Livinestone et al. (2006) near MJD 52210 ± 10 

^ Phase coherence lost over 78 d data gap starting at MJD 52837 (see also Livingstone et al.. 2006) 

^ Phase coherence lost after MJD 53880 (May 24, 2006) during the magnetar-like outburst (Gavriil et al.. 2008) 



4. Timing 

4. 1 . Timing solutions: ephemerides 

RXTE observed PSR J 1846-0258 for the first time in the period 
April 18-21, 1999 (MJD 51286-51290) for about 39 ks. The ob- 
servation was split in 1 1 short observations (sub-observations) 
of duration ranging from 2.5 to 19.9 ks. We barycentered the 
PCA event arrival times using the Chandra X-ray observatory 
(CXO) sub-arcsecond position of PSR 11846- 0258, {a,S) ^ 
a8H6"24^.94,-02°58'30'.'l) for epoch J2000 dHelfand et al.L 
l2003h . which con-esponds to (l,b)=(29.712015,-0.240245) in 
Galactic coordinates. In each of these sub-observations a co- 
herent pulsed signa l at a rate of ~ 3.0827 Hz (see also 
iGotthelf et al.L l2000l) could be detected in the bary centered time 
series. We used a Z^-test (Buccheri et al.L 119831) search in a 
small, typically 5 independent Fourier steps (Avips = 1/t, 
in which t represents the time span of the data period), win- 
dow around the predicted pulse frequency. The restricted search 
yielded for each sub-observation a best estimate of the rotation 
rate at the gravity point of the sub-observation. Because the Z^- 
test is distributed as a for 2x1 degrees of freedom a Icr 
error estimate on this optimum value can easily be derived by 
determining the intersection points of the measured Z^-test dis- 
tribution near the optimum with the value of Z^^^^^^ - 2.296. The 
set of optimum pulse frequencies and Icr uncertainties versus 
time can be used to obtain an incoherent timing solution v,„c(f) 
for the rotation behaviour of the pulsar Phase folding the arrival 
times of each sub-observation on the appropriate optimum fre- 
quency results in pulse profiles for each sub-observation. In or- 
der to derive the mutual phase relation of the pulse maxima we 
use cross-correlation techniques to obtain a high-statistics pulse- 
profile template. This template, based on 1 1 sub-observations, 
was subsequently used in the time of arrival (TOA) determina- 
tion procedure, described below. 



The accuracy of the timing parameter^ (v, v, v) can be im- 
proved considerably relative to the incoherent set by maintaining 
pulse-phase coherence over the entire time interval of the con- 
sidered set of sub-observations. This requires an accurate de- 
termination of pulse arrival times for each sub-observation in 
the time interval considered. These times are obtained by cross- 
correlating the pulse profiles of each sub-observation with the 
high-statistics template. The global maximum in the correlation 
diagram determines the time shift (AOmax/i') to be applied to 
the chosen time zero point of the sub-observation (normally the 
gravity point) to align it with the template. The global maxi- 
mum of the correlation function is derived by fitting a truncated 
Fourier series (typically 5 harmonics are used) to the measured 
correlation function, evaluated at discrete step points, in order to 
suppress fluctuations. The latter could be significant in case of 
weak pulsed signals. An error estimate 5<t) on AO^ax has been 
obtained by bootstrapping both the pulse-profile of the active 
sub-observation and that of the template assuming an underly- 
ing Gaussian distribution. Typically 1000 resamplings are used 
in subsequent correlation analysis. The Icr error on AO^ax corre- 
sponds to the width of the distribution of the correlation maxima. 
Thus, we obtained a set of n pulse-arrival times with correspond- 
ing error estimates, (fJJ, Afp). Next, the arrival times were folded 
upon the incoherent or any other reasonable timing solution to 
obtain the residuals 

A<D^ = Vine ■ (fp - fp) + \vi„c ■ (4 - fp)' + ■ (4 - 0' 

Finally, we determined the increments Av, Av and Ai> which 
minimize \h&x^ statistics 2t=i (1^)^ taking into account the sta- 
tistical weights of the individual points. This approach yields a 
phase-coherent (or phase-connected) timing solution (Vc, Vc, Vc) 

' In this work we limited ourselves to a maximum of three timing 
parameters or less in case of small time spans. 
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describing every revolution of the pulsar over the time span for 
which the solution is valid. 

As of January 30, 2000 RXTE monitors PSR J 1846-0258 
regularly, and from observations performed till November 7, 
2007 (MJD 51574-5441 1) we derived phase coherent timing so- 
lutions (ephemerides) when possible. Table |2] lists the current 
set of phase-coherent ephemerides, as derived in this work us- ' 
ing the method described above, including the one (entry # 0) ^ 
covering the four days of the very first RXTE observations of 
PS RJ1846-0258. 

[Livingstone et alj (l2006l) presented equivalent timing solu- 
tions for the range MJD 51574-53578, covering about 5.5 year 
of X-ray timing data. Our timing results cross the magnetar-like 
radiative outb urst event that occurred between May 24- 3 1 , 2006 
(see e.g. Gavriil et al., 2008 HKumar & Safi -Hai-b, 2008). Entries 
1 till 5 of Table |2] cover the so-called pre-outburst period, which 
is characterized by very stable and highly predictable timing be- 
haviour allowing an accurate measurement of the breaking in- 
dex of the pulsar spin-down of 2.65 + 0.01 ( [Livingstone et al.L 
12006 ). Our timing results for the pre-outburst period (MJD 
51574-53880) ar e fully compatible with those obtained by 
[Livingstone et al.l ([20061) for a slightly smaller time span (MJD 
51574-53578). It demonstrates that PSR J1846-0258 continued 
to slow-down very predictably from MJD 53578 till MJD 53880, 
the date of the last sub-observation before the magnetar-like out- 
burst. 
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Fig. 2. The rotation behaviour of PSR J 1846-0258 from Jan. 30, 
2000 to Nov. 7, 2007 as determined from RXTE monitoring ob- 
servations relative to the pre-outburst timing model (entry # 5 
of Table |2]i. Note the dramatic change during the radiative out- 
burst and its subsequent decay (e-folding time scale ~ 55 days). 
The presence of a huge spin-up glitch near MJD 53883 + 3 is 
clearly visible. Phase coherence is caught-up again as of MJD 
54126. Lines: phase-coherent timing models. Data points with 1- 
cr errors: incoherent solutions. 



4.2. A small and a major spin-up glitch 

We also confirm th e loss of phase-coheren ce beyond MJD 
52237, attributed by [Livingstone et akl ([20061) to a small glitch 
of size Av/v = 2.5(2) x 10"^ and Av/v ~ 9.3(1) x lO""* near 
MJD 52210 + 10. Note, however, that the TOA for the sub- 
observation on MJD 52237 can still perfectly be predicted by 
the timing parameters of entry # 1 of Table |2] indicating that the 
small glitch very likely occured beyond MJD 52237, but before 
MJD 52253, the next sub-observ ation, thus later than the glitch 
epoch of 52210 + 10 reported bv [Livingstone et all ([2006h . 

As of the onset, beyond MJD 53880, of the radiative out- 
burst till MJD 54126 (January 26, 2007) phase coherence is 
completely lost du e to highly increased timing noise (see also 
[Gavriil etali[200a ). The sparse sampling in this period and the 
high timing noise exclude the construction of a phase coherent 
solution during the outburst and its subsequent decay except for 
a small period of 55 days from MJD 53978 till 54033 (entry #6 
of Table|2]). Phase coherence is caught-up again on MJD 54126 
(entries #'s 7 and 8 of Table |2|. Note the very high value of v 
and the increased |v| values of entries #'s 6 and 7. All timing in- 
formation both coherent and incoherent, if the former is lacking, 
is shown in Fig.|2] In this figure the frequencies (coherent, solid 
lines; incoherent, data points) are relative to the timing model 
valid just before the radiative outburst (entry # 5 of Table|2]i. 

From Fig.|2]it is clear that the radiative outburst was accom- 
panied by a huge spin-up glitch near MJD 53883 + 3 as sus- 
pected already bv lGavriil et alj ([2008 h, but now clearly demon- 
strated. Depending on the number of sub-observations used in 
the construction of an incoherent timing model since the on- 
set of the outburst we obtain a glitch size Av/v in the range 
(2.0 - 4.4) X 10"^. This fractional-frequency-jump size ranks in 
the top of the glitch-size distribution for rotation-powered pul- 
sars, but is typical for AXP glitches (Dib et al., 2008), although 
the number of observed AXP glitches is still too small to se- 
curely establish the underlying distribution. 
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Fig. 3. The pre-outburst (Revs. 49-441) pulse profile of 
PSR J1846-0258 as observed by IBIS ISGRI in the 20-150 keV 
energy band with a detection significance of 9.6cr. Error bars 
are Icr on measured counts. Superposed is the best fit model 
composed of a PCA-template profile for the ~ 2.9-8.3 keV band 
(dashed line) and a DC-level (long dashes; horizontal line) with 
its Icr-error estimates (see Sect. 15. lb . There is no significant dif- 
ference in shape. 



4.3. Pre-outburst pulse profiles: INTEGRAL ISGRI, RXTE 
PCA and HEXTE 

The timing analysis of INTEGRAL ISGRI data consists of se- 
lecting, screening and finally phase folding barycentered time 
series with appropriate ephemerides. The observation and event 
selection procedures are identical to those applied for the imag- 
ing studies of ISGRI data. In order to reduce in the timing 



L. Kuiper and W. Hermsen: High-energy characteristics of the schizophrenic pulsar PSR J1846-0258 in Kes 75 



7 




0.0 0.5 1.0 1.5 2.0 0.5 1.0 1.5 2.0 0.5 1.0 1.5 2.0 

Pulse Phase 



Fig. 4. Collage of pulse profile of PSR J1846-0258 using pre-outburst data from RXTE PCA (3-35 keV: a-c), HEXTE (15-250 
keV: d-f) and INTEGRAL ISGRI (20-300 keV: g-i). The pulse signal is detected for the first time up to ~ 150 keV making 
PSR J1846-0258 the fifth pulsar showing up at soft y-rays. The pulse profile shape is energy independent: a broad single-peak 
asymmetric pulse. Error bars are Icr on measured counts. 

analysis the backgound level consisting of non-source counts, 
only detector pixels are considered which have been illumi- 
nated by the source for at least 25%. In this work, we further 
screen the data on short duration (lasting less than a second) 
bursts, particularly those originating from the soft gamma-ray 
repeater SGR 1806-20, which pollute the genuine pulsed sig- 
nal from PSR J1846-0258. Then, the on-board-registered event- 
time stamps of the selected events are corrected for known in- 
strumental (fixed), ground station and general time delays in the 
on-board -time versus Terrestrial-Time (TT or TDT) correlation 
(see e.g. IWalter et all 12003 ). The resulting event times in TT 
are barycentered (using the JPL DE200 solar system ephemeris) 
adopting the CXO position of PSR J 1846 -0258 and the instanta- 
neous INTEGRAL orbit information. These barycentered events 
(in TDB time scale) are finally folded upon the appropriate tim- 
ing model composed of v, v, v and the epoch fo, as obtained in the 
phase-coherence analysis of RXTE PCA data (see Table|2]i. The 
TDB-time to pulse-phase conversion taking into account con- 
sistent profile alignment by subtracting (Sq is provided by the 
following formula: 

1 ,1 , 

<D(/) = V ■ (f - fo) + -V ■ (f - fo) +-v-(t- toY - 3)0 
/ 6 

Thus, we produced pulse-phase distributions for differential en- 
ergy bands of width 1 keV between 15 and 300 keV. For 
the "stable-timing-behaviour" period before the onset of the 
magnetar-like outburst we obtained a 9.6cr signal (Zj-test) in the 
20-150 keV band (see Fig. [3]) combining all timing data from 
INTEGRAL observations performed between revolutions 49 and 
441 (see Table[T]l. The total Good-Time-Interval (GTI) exposure 
of the 2599 Scw's included in this combination amounts ~ 6 Ms, 
which translates in an effective on-source exposure of ~ 3 Ms. 



The right hand panels g, h and i of Fig.|4]show the ISGRI pulse 
profiles in the three differential energy bands, 20-35 keV, 35-60 
keV and 60-150 keV. The pulsed-signal detection significances 
are 5.6cr, 7.0cr and 4.2cr, respectively, making PSR J1846-0258 
the fifth pulsar detected in the soft y-ray band, after the Crab, 
Vela, PSR B 1509-58 and PSR B0540-69 pulsars. 

We wish to compare the ISGRI pulse-phase distributions 
with the time-averaged HEXTE and PCA pulse profiles from 
all RXTE observations performed before the magnetar-like out- 
burst. The timing analysis of HEXTE data, covering an energy 
window c omparable to ISGRI , is equivalent to the approach fol- 
lowed bv iKuiper et all (l2006h (see Sect. 3.1 of that paper). We 
collected 344.122 ks and 335.770 ks exposure, both corrected 
for dead time and reduction in efficiency due to off-axis obser- 
vations, for HEXTE cluster A and B, respectively. Phase folding 
HEXTE barycentered-time series with the ephemerides derived 
from the simultaneously taken PCA data (see Table |2]i yields 
the high-energy phase distributions from HEXTE. Differential 
HEXTE profiles are shown in panels d-f of Fig.fflforthe 15.6-28 
keV, 33.1-60.1 keV and > 60.1 keV energy bandfl respectively. 
The corresponding significances are 1 l.Ocr, 8.7cr and 4.3cr. The 
ISGRI and HEXTE profile shapes are fully consistent. Finally, 
we also included in Fig. |4] the PCA pulse profiles from pre- 
outburst observations for energies between ~ 2.9 and ~ 32.5 
keV. As can be seen, the morphology of the pulse profile for 
energies above ~ 2.9 keV appears to be energy independent: 
a broad single-peak asymmetric pulse with a somewhat steeper 
rise than fall. 



- Spectral data from the 28-33.1 keV HEXTE band have been ig- 
nored because of the presence of a huge background line due to the 
activation of iodine. 
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Fig. 5. Comparison of the pulse-profiles accumulated in the PCA 
PHA band 7-19 (~ 2.9 - 8.3 keV) during the early phase of the 
outburst (Outburst A1H-A2, see Table |3] solid histogram with 
Icr error bars) and from the pre-outburst observations (high- 
statistics dotted histogram). No significant shape change is ob- 
served, see text. The Y-axis specifies the number of counts per 
bin. 

4.4. Pulse-profile morphology during the outburst 

To investigate whether the radiative outburst, most pronounced 
at soft X-rays, was accompanied by a pulse-morphology change 
we compared the high-statistics pre-outburst PCA profile with 
the profile from RXTE-PCA observations during the early phase 
of the outburst. For the latter we used data collected between 
MJD 53894 and 53970 (Outburst Al + A2, see Tableland also 
Fig. Q, omitting those observations in which the reported five 
short (0.1 s) magnetar-like bursts occurred. We selected only 
events from PHA channels 7 to 19, roughly corresponding to 
2.9-8.3 keV. Because phase coherence was lost after the onset 
of the outburst, we generated the "outburst" profile by stacking 
the properly shifted profiles from individual observations apply- 
ing a cross-correlation analysis which yields the required phase 
shifts (see eg. Section 3.1 of Kuiper et al., 2004, for a similar 
analysis). Fig. |5] shows the superposed aligned pre-outburst and 
outburst profiles. We fitted the outburst profile in terms of a con- 
stant and the shape of the pre-outburst profile with a free scale. 
The absolute values and signs of the deviations from the best fit 
were used in a combination of two independent statistical tests, 
a Pearson test and a run test. The combined probability of 
these tests is 3.2%, i.e. both profiles are different at a 2.15cr sig- 
nificance level, indicating that there is no significant change in 
shape between the pre-outburst and outburst profile. 

5. Pulsed Emission Spectra 

5. 1 . Analysis methods 

From the pulse-phase distributions, F{(t>, Ep^j^), derived for 
RXTE PCA, HEXTE and INTEGRAL IBIS ISGRI, we ex- 
tracted pulsed excess counts by fitting two different model func- 



tions to the measured pulse-phase distributions in selected en- 
ergy bands: 

1) a truncated Fourier series with - 1 harmonics, 
A?((D) = fl() + Ok cos(2;ryt(D) + sin(2;rfc<D) 

2) a PCA-template model, A^(<D) ^a + bx r((D) 

For model 1, a number of = 3 was sufficient to adequately 
describe the measured profiles. The template function T{(S>) used 
for model 2 is based on a high-statistics PCA pulse profile for 
energies between ~ 2.9 - 8.3 keV. 

The fit function minimum plus its error are subsequently 
used to determine the number of pulsed excess counts i.e. the 
number of counts above this minimum (=unpulsed or DC level) 
along with an error estimate. For pulse profiles with a strong 
pulsed signal both methods yielded consistent results, however, 
for weak pulsed signals method 1 overestimates the number of 
pulsed excess counts slightly (a small positive bias). Therefore, 
we have used the template fit model subsequently, assuming thus 
that the genuine underlying pulse profile did not vary with en- 
ergy which is a very reasonable approximation (see e.g. Fig.[3]l. 

For the PCA we constructed time-averaged energy response 
matrices for each PCU separately taking into account the differ- 
ent (screened) exposure times of the involved PCU's during the 
time period of interest. For this purpose we used the f tools ver- 
sion 6.4 programs pcarsp and addrmf. To convert PHA channels 
to measured energy values, fipHA, for PCU combined/stacked 
products we also generated a weighted PCU-combined energy 
response matrix. 

For HEXTE we employed cluster A and B energy-response 
matrices separately, taking into account the different screened 
on-source exposure times and the reduction in efficiency in case 
of off-axis observations. The on-source exposure times for both 
clusters have been corrected for considerable dead-time effects. 

We assume simple power-law models, Fy - K ■ {EylE^y^ 
with r the photon-index and K the normalization in ph/cm^s keV 
at the pivot enemy Eq, for the underlying photon spectra. In case 
of RXTE PCAQand HEXTE these models have been fitted in a 
forward folding procedure using appropriate response matrices 
to obtain the optimum spectral parameters, K and F, and the re- 
constructed spectral flux points from the observed pulsed count 
rates. 

A different approach has been used for IBIS ISGRI: For each 
energy band we scaled the derived pulsed excess counts by the 
number of pulsed excess counts extracted for the Crab pulsar in 
exactly the same measured energy window. These ratios were 
subsequently multiplied by the pulsed photon spectrum of the 
Crab as derived from deep HEXTE observations (see Sect. 3.4 of 
iKuiper et al.L l2006l for details) taking into account the different 
effective on-source exposures for both PSR J 1846-0258 and the 
Crab. 

5.2. Pulsed pre-outburst 3-300 keV spectrum 

The reconstructed pulsed emission spectra from pre-outburst 
RXTE and INTEGRAL observations are shown in Fig.|6]as aqua 
coloured data points for the PCA (~3-30 keV), blue data points 
for HEXTE (15-250 keV) and magenta coloured symbols for 
INTEGRAL ISGRI (20-300 keV). The individual PCA, HEXTE 

^ An absorbing Hydrogen column of Nh = 3.96x 10^^ cm~^ h as been 
assumed for the PCA spectral analysis fsee lHelfand et alil2003h and in 
all PCA spectra shown in this paper the interstellar absorption has been 
modeled out. 



L. Kuiper and W. Hermsen: High-energy characteristics of the schizophrenic pulsar PSR J1846-0258 in Kes 75 



9 



10-": 



T 



1 0" 




Pulsed (INTEGRAL IBIS-ISGRI) 
Total (INTEGRAL IBlS-lSGRl) 
PSR (Chandra ACIS) 
PWN (Chandra ACIS) 
Pulsar + PWN (Chandra ACIS) 
Pulsed (RXTE HEXTE) 
Pulsed (RXTE PCA) 



10-n 



K [McV] 

Fig. 6. Total and pulsed (unabsorbed) high-energy spectrum (~ 
1 - 300 keV) of Kes75 / PSR J 1846-0258 from pre-outburst 
RXTE PCA and HEXTE (both pulsed) and INTEGRAL IBIS 
ISGRI (pulsed and total) observations in an E^F representa- 
tion. The spectral coverage is extended towards 1 keV using the 
Chandra spec t ral-fit results (1-7 keV; unabsorbed) obtained by 
iHelfand et all (|2003|) for the total emission spectrum of the pul- 
sar and the PWN. Note that both the PWN contribution and DC 
emission from PSR J1846-0258 drop to undetectable levels at 
energies near 150 keV: the total emission is consistent with be- 
ing 100% pulsed above ~ 150 keV. Error bars on data points are 
lo-. 

and ISGRI spectra are mutually consistent in overlapping energy 
bands. 

A single power-law model fit to the PCA data points (~ 3-30 
keV) alone gave an optimum photon index F of 1.161 + 0.004. 
However, the obtained reduced of 2.38 for 13 degrees of 
freedom is poor, and fitting a power-law with an energy de- 
pendent index, a so-called "curved power-law model" Fy - 

K • (iiy/Eo) provided a significant improvement of 

3.4cr over the simple power-law model, mainly due to better de- 
scribing the flux measurements below 4.5 keV. 

In order to obtain an accurate description of the pulsed spec- 
trum over two decades in energy we also fitted a power-law and 
curved power-law model over the full ~ 3-300 keV energy band 
to the combined PCA, HEXTE and ISGRI pulsed flux measure- 
ments. The best-fit power-law model, shown as a solid black line 
inFig.|6l has aphoton indexFof 1.20 + 0.01, slightly softer than 
the fit to just the fluxes below 30 keV. The energy fluxes of the 
pulsed emission (unabsorbed) in the 2-10 keV and 20-100 keV 
bands are (2.38 + 0.02)x lO^'^ erg/cm^s and (15.2±0.14)x lO^^^ 
erg/cm^s, respectively. Assuming a 1 steradian beam size for the 
pulsed emission and a distance of 10 kpc the luminosities in the 
2-10 and 20-100 keV bands are L^-io = 2.27(2) x 10^^ ■ rj ■ 

erg/s and Lf-^^ = 1.45(1) x 10^"* ■ '7 ■ -^^o ^^g/s (// is the beam- 
size), respectively. These numbers convert into the following 
efliciencies (Lx/E.d), 4"'° = 0.027% ■ 77 ■ 4, and 4°-'°° = 
0.18% • J] ■ respectively. The fit with a curved power-law 
model was again better (4.9cr improvement over the power-law 



model; it is shown in Fig.fTTTi for the following best-fit parameter 
values: K = (8.40 + 0.12)x lO^^'ph/cm^s keV, f = 1.295 + 0.016 
and 6 = 0.096 + 0.014 for a pivot energy of 17.6209 keV. 



5.3. Pulsed fractions 

In Fig. |6] the total (i.e. PSR J1846-0258 pulsed plus un- 
pulsed/DC and PWN contribution) pre-outburst ISGRI spectrum 
(see Sect. 13.21 ) is presented as purple data points across the 20- 
300 keV band as well as the best fit power-law model to these 
measurements (purple dashed line). Using the energy fluxes of 
pulsed and total emission in the 20-100 keV band we could de- 
rive a lower-limit on the pulsed fraction of 44% in the 20-100 
keV band. Moreover, it is evident that near 150 keV the pulsed 
emission becomes consistent with the total emission i.e. 100% 
pulsation. As a result, both, the spectra of the DC-emission from 
the pulsar and of the PWN must bend down between 20 and 100 
keV to significantly lower levels. 

To extend the spectral coverage to lower energies we in- 
cluded in Fig. |6]the spectral model fits obtained from Chandra 
data in the soft X-ray ba nd (1-7 keV) by Helfand et al. (200|) 
(see for recent analyses iKumar & Safi-Harb. ,2008; .Ng et afl. 
I2OO8I) for the total emission from PSR J1846-0258 (F = 1.39 + 
0.04; pulsed plus unpulsed/DC; solid orange line) and the PWN 
(F - 1.92 + 0.04; dashed orange line), separately. We also show 
the sum of these two models as a dotted orange line represent- 
ing the total emission from PSR J1846-0258/PWN in the 1- 
7 keV band, to be compared with the total emission as mea- 
sured by ISGRI above 20 keV. Dividing the pulsed emission es- 
timated in the 0.5-10 keV band through extrapolation to lower 
energies of the ~ 3-300 keV pulsed spectrum, (2.98 + 0.03) x 
10 erg/cm^s, b y the total poi nt-source emission 9.50 x 10"'^ 
erg/cm^s (HelfandHifl l2003h from PSR J 1846-0258 in the 
Chandra 0.5-10 keV band we obtained an average pulsed frac- 
tion of 31%. Therefore, comparing this value with the lower- 
limit of 44% derived for the 20-100 keV band, the pulsed frac- 
tion is steadily growing with energy to become within the large 
errors consistent with ~ 100% near 150 keV. This is primarily 
an indication that the DC spectra have to bend down. 



5.4. Spectral evolution (3-30 keV) before, during and after 
the outburst 

The X-ray spectral index of the total emission from 
PSR J 1846-0258 (pulsed plus DC emission) measured with 
Chandra between 1 a nd 7 keV ch anged from a pre-burst value 
in 2000 of F ~ 1.35 jHelfand et al., 2003; Kumar & Safi-HarR 
20 08D to a value of F ~ 1.93 (IKumar & Safi-Harbl I206i 
Gavriil et al 112Q08) during the magnetar-like outburst in 2006. 
In order to better characterize the spectral evolution we took ad- 
vantage of the high-statistics 3-30 keV PCA data to derive spec- 
tra for the pulsed emission before, during and after the outburst. 

Fig. I2] indicates in the top panel the selected time seg- 
ments for the PCA analysis in comparison with the count- 
rate lightcur ve (2-60 keV) of PSR J 1846-0258 published by 
iGavriil et al.l (j2008): Segment Pre-II covers an interval just be- 
fore the outburst; segments Al, A2 and B cover the reported out- 
burst; Post-I to Post-Ill are time intervals after the outburst for 
which the count-rate lightcurve indicates that PSR J 1846-0258 
returned to its pre-burst flux level. Table [3] shows the details of 
the PCA time segments including the identifiers of the first and 
last sub-observation along with the corresponding start and end 
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Table 3. Time intervals used in the spectral analysis of the RXTE PCA data. 
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Obs. id. 


MJD 

Begin/End 





Screened PCU exposure [ks] 
1 2 3 


4 


Name 


40140-01-01-01 
90071-01-10-00 


51286 
53340 


940.624 


217.960 


940.608 


814.520 


350.200 


Pre-outburst I 


90071-01-12-00 
92012-01-12-00 


53395 
53879 


203.712 


33.768 


207.944 


107.176 


13.120 


Pre-outburst II 


92012-01-14-00 
92012-01-19-00 


53894 
53926 


20.064 


4.880 


24.008 


13.712 


2.832 


Outburst Al' 


92012-01-20-00 
92012-01-25-00 


53935 
53970 


11.096 


9.616 


20.896 


6.720 


2.992 


Outburst Ai^ 


92012-01-26-00 
92012-01-32-00 


53978 
54019 


23.032 


8.232 


31.176 


10.176 


9.664 


Outburst B 


90071-01-18-00 
91071-01-16-00 


54126 
54160 


32.096 


3.952 


39.632 


7.568 


13.064 


Post-outburst I 


91071-01-17-00 
93010-01-05-00 


54165 
54313 


78.048 


11.330 


91.360 


22.192 


32.440 


Post-outburst II 


93010-01-06-00 
93010-01-19-00 


54319 
54411 


52.688 


12.496 


72.888 


18.392 


18.792 


Post-outburst III 



'Obs.id. 92012-01- 
^Obs.id. 92012-01- 



13-00 on MJD 53886 has been excluded because of the occurence of 4 SGR-like bursts 
21-00 on MJD 53943 has been excluded because of the occurence of 1 SGR-like burst 
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Fig. 8. Evolution of the (unabsorbed) pulsed 
spectrum (-3-30 keV) of PSR J1846-0258 as 
derived from RXTE PCA data. Note the dra- 
matic softening of the spectrum during the 
early phase of the outburst (purple data points) 
and the gradual return to its pre-outburst spec- 
trum (black data points). Also indicated for 
each time segment are the power-law model fits 
obtained from > 9 keV data. Data points are 
given with Icr errors. 



time in MJD. Also, the screened exposure time per PCU is listed 
together with the name of the time segment. 

The bottom panel of Fig. |7] shows that we derived for each 
of these time segments coherent or incoherent timing solutions, 
allowing us to produce pulse profiles and subsequently spectra of 
the pulsed emission for each segment. Comparing the upper and 
bottom panel of Fig. |2l it is clear that the major glitch triggered 
the radiative outburst. 

For each PCA observation-time segment we had sufficient 
counting statistics to produce pulsed spectra (~3-30 keV). Fig. 
[8] shows the (unabsorbed) pulsed spectra for all time segments 
listed in Table |3] including the multi-year-average spectrum of 
Pre-outburst I, but combining the three post-outburst spectra. 
The latter three spectra are combined, because these are statisti- 



cally identical and the individual flux values could not be distin- 
guished when plotted in this figure. It is immediately evident that 
the pulsed spectrum evolved smoothly in flux and shape during 
the burst. A dramatically softer pulsed spectrum was measured 
directly after the glitch (time interval Al, purple data points). 
Obviously, a strong soft component was added to the pre-burst 
hard (F ~ 1.16) non-thermal spectrum, totally dominating the 
pulsed emission below ~ 10 keV. During time segment Al the 
pulsed (energy) flux for energies 2-10 keV increased by a fac- 
tor ~ 5. This soft component quickly reduced in intensity during 
the outburst (A2, red; B, orange data points), and is not detected 
anymore in the post-outburst spectra which are within statistics 
identical to the Pre-outburst I and II spectra. In Fig. |8] one can 
see that the time variability is less pronounced for energies above 
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Fig. 7. Pulsed-flux lightcurve (P CA 2-60 keV) arou nd the out- 
burst (top panel) adapted from iGavriil et alj (l2008h . The pre- 
outburst averaged pulsed-flux level is shown by the horizontal 
dashed line along with its tiny Icr error (plotted at MJD 53400). 
The coverage of the INTEGRAL observations near the outburst 
period are indicated with revolution-number intervals in the up- 
per part of the top panel. In the lower part of the top panel are 
shown the time segments used in the spectral-evolution study 
of the pulsed signal using 3-30 keV PCA data. The lower panel 
displays the pulse frequency evolution across the outburst period 
(zoom-in of Fig.|2|. The phase-coherent solution in the time in- 
terval of 55 days around MJD 54000 is indicated. Errors are Icr. 
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Fig. 9. Evolution of the pulsed flux of PSR J 1846-0258 in the 
10-30 keV band as derived from a power-law fit to RXTE PCA 
> 9 keV spectra (photon index fixed to its pre-outburst value 
of 1.16). Error bars are Icr. The average flux level in the Pre- 
outburst 1 time segment (duration ~ 5.6 yrs) is indicated with the 
broken line together with its +lcr uncertainty (dotted fines). 

have larger statistical errors and can be fitted with single power- 
laws. However, when fitted with two power-laws, these are also 
consistent with the sum of two spectral components, a soft one 
with index F, = 3.27 + 0.27 (segment A2) and = 4.74 + 0.71 
(segment B), along with F/, = 1.16. Note the apparent gradual 
softening of the soft component below ~ 10 keV during the de- 
cay of the outburst, indicative for cooling. 

In Fig. [To] we present the variation of the measured pulsed 
fluxes below 10 keV. These are derived using the above two- 
power-law fits in the outburst time segments, and the single 
power-law fit (index 1.16) for pre- and post-outburst. The decay 
of the excess photon flux measured above the non-thermal pre- 
outburst flux level (indicated in Fig. [TOl i is consistent with the 
ex ponential decay with 1/e time constant of ~ 55.5 day reported 
bv lGavriiletan(l2008h . 



~ 9 ke V. In fact, we verified that for all selected time intervals the 
spectra above 9 keV are fully compatible with a single power- 
law model shape with the pre-outburst index F - 1.16. In Fig. 
[8] these power-law model fits, all with a fixed photon index F of 
1.16, are superposed for each time segment. The corresponding 
pulsed fluxes (10-30 keV) are shown in Fig.|9l including also the 
individual data points for the three post-outburst time segments. 
We see, that also in this case the flux enhancement was maximal 
directly after the glitch in time segment Al, but amounts only ~ 
35%, which represents a 4.3cr enhancement. 

It is evident that during the outburst the spectral shape over 
the 3-30 keV band varied. A single power-law model fit for time 
segment Al gave a poor fit (F = 2.46 + 0.01 with a reduced 
X^r - 1-69 for 13 degrees of freedom). A model fit with two 
(free) power-laws rendered a 3.2cr improvement. An excellent 
fit is obtained with a soft power-law index Fj = 2.96 + 0.06 
along with a hard photon index of F/, = 1.16, fixed at the pre- 
outburst value. The pulsed spectra in time segments A2 and B 



6. Summary 

In this paper we have presented detailed high-energy charac- 
teristics of the enigmatic PSR J 1846-0258 and its PWN, us- 
ing the complementary spatial, spectral and timing capabili- 
ties of RXTE PCA and HEXTE and INTEGRAL'S IBIS/ISGRI. 
Particularly important is to see how the characteristics com- 
pare before, during and after the unique magnetar-like out- 
burst in order to recognize constraints which are important for 
the interpretation of the transient magnetar-like behaviour of 
PSRJ1846-0258. 

6.1. Total spectra above 20 keV of PSRJ1846-0258/Kes 75, 
before, during and after outburst 

1 - Using all available INTEGRAL data taken with ISGRI 
during 2003-2006 before the magnetar-like outburst in June 
2006, adding up to 3.0 Ms effective on-source exposure, the 
time-averaged 20-300 keV spectrum for the total emission of 
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Fig. 10. Evolution of the pulsed flux of PSR J 1846-0258 in the 
2-10 keV band as derived from RXTE PCA data. Error bars are 
Icr. The average non-thermal flux level in the Pre-outburst 1 time 
segment (duration ~ 5.6 yrs) is indicated with its ±lcr uncer- 
tainty like in Fig.|9] but the levels practically coincide. 

PSRJ1846-0258/Kes 75 can be represented by a power-law 
with photon index F = 1.80 + 0.06 and an energy flux (20-300 
keV) of (6.62 + 0.35) x IQ-" erg/cm^s {Sect.\3^. 

2 - ISGRl observed PSR J1846-0258/Kes 75 unfortunately only 
towards the end of the outburst for 214 ks. No change in spectral 
shape was observed (F - 1.75^^ 3 j), but there is an indication 
at the 2.3cr level for an enhanced flux level (by ~ 52%) during 
outburst decay (Sect. [Ob . 

3 - About one year after the onset of the outburst, the measured 
(415 ks) total post-outburst spectrum in the 20-300 keV band is 
within Icr errors of ~ 15% consistent in shape and flux with the 
pre-outburst spectrum (Sect. \3.4[ . 

6.2. Ephemerides of PSR J1 846-0258 before, during arid 
after outburst; glitches 

1 - Exploiting the multi-year monitoring observations of RXTE 
and the high count rate of the PCA we derived phase-coherent 
ephemerides for the 6.5 years be fore the outburst, consis- 
tent with the results reported by iLivingstone et al.l (l2006h . 
PSR J 1846-0258 was manifesting itself as a stable young high 
B-field pulsar with a breaking index of the spin-down of 2.65 + 
0.01 (Sect.\4j\ TaWelll. 

2 - During the ou tburst phase coherence was lost (see also 
iGavriil et al.L '2008), except for a short period of 55 days from 
MJD 53978 till 54033 in which we found a phase-coherent solu- 
tion. Nevertheless, incoherent solutions have been derived over 
the full duration of the outburst. Phase coherence was caught-up 
again after the burst on MJD 54126, when the PCA (2-60 keV) 
pulsed flux was back to its average pre-burst value (Sect. 14.21 
TabJeH Fig.\3. 

3 - We showed that the onset of the radiative outburst was 
accompanied by a huge spin-up glitch near MJD 53883 + 3 
with a glitch size Av/v in the range (2.0 - 4.4) x 10"^. This 
fractional-frequency-jump size ranks in the top of the glitch- 
size distribution for rotation-powered pulsars. Furthermore, we 



confirm the occurrence o f a small glitch near MJD 52210 + 10 
(iLivingstone et al.Ll2006 l), however, for a slightly different epoch 
somewhere between MJD 52237 and 52253 (Sect.^ 

6.3. Pulse profiles 

1 - Using the high statistics of all pre-outburst observations, 
ISGRI and RXTE HEXTE both measure significant and con- 
sistent pulse profiles up to 150 keV. Namely, a broad single 
asymmetric pulse, making PSR J 1846-0258 the fifth rotation- 
powered pulsar detected at hard X-rays/soft gamma rays (Sect. 
[QFig.Dll. 

2 - Comparison of the pulse profiles measured with the PCA, 
HEXTE and ISGRl between 2.9 and 150 keV shows that the 
pulse shape does not change with energy (Sect. 14.31 Figs. [3] & 
13. 

3 - The shape of the profile measured with the PCA between 
2.9 and 8.3 keV during the magnetar-like outburst, when the soft 
component dominated the emission, exhibits a shape which is 
fully consistent with the high-statistics pre-outburst pulse profile 
in the same energy band. This is evidence for a stable geometry 
(Sect. m Fig. 

6.4. Pre-outburst pulsed spectra and pulsed fractions; 
PWN-DC spectrum 

1 - The time-averaged spectra of the pre-outburst pulsed emis- 
sion as measured with the PCA, HEXTE and ISGRI between 
~ 3 keV and 300 keV are consistent in overlapping energy 
bands. The best-fit power-law model has a photon index F = 
1.20 + 0.01. The energy fluxes of the pulsed emission in the 2-10 
keV and 20-100 keV bands are (2.38 + 0.02) x lO^'^ grg/cm^s 
and (15.2 + 0.14) x 10"'^ erg/cm^s, respectively. Fitting a curved 
power-law model to just the PCA spectrum (3-30 keV) and to 
the total combined pulsed spectrum (3-300 keV) gives a signifi- 
cantly better fit. (Sect.\5i^ Figs.[6l&f771l. 

2 - Comparing the measured pulsed flux (only PSR J 1846-0258 
pulsed) and the total (PSR J 1846-0258 pulsed + DC + PWN) 
measured by ISGRl in the band 20-100 keV gives a lower limit 
to the pulsed fraction of 44% for the emission from the compact 
object. Around 150 keV the measured pulsed and total fluxes are 
consistent with a pulsed fraction of 100%, however the statisti- 
cal uncertainties are large at these energies. Comparing the PCA 
pulsed flux with the total point-source (PSR J 1846-0258 pulsed 
+ DC) flux measured with Chandra, gives an estimated pulsed 
fraction of 31% for the 0.5-10 keV band (Sect.\53i Fig.\^). 

3 - The above discussed variation of pulsed fraction with energy 
implies that the high-energy spectrum of the PWN, as well as 
that of the point-source DC emission of PSR J1846-0258 mea- 
sured by Chandra below 7 keV will strongly bend down in the 
INTEGRAL energy window above 20 keV to drastically lower 
(£¥) levels at 100 keV (Sect.O Fig.\^. 

6.5. Spectral evolution (3-30 keV) during magnetar-like 
outburst 

1 - In the first 32 days after the major glitch (time segment 
Al) the 3-30 keV pulsed spectrum can be represented with two 
power-law models, a soft component with index F, - 2.96 + 0.06 
and a hard component with the pre-outburst value F/, - 1.16. The 
2-10 keV (energy) flux increased by a factor ~ 5 and the 10-30 
keV (energy) flux increased with only 35% (Sect. \5~4\ Figs.^& 
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2 - Above 9 keV all spectra during outburst and also the 3-30 
keV pre- and post-outburst spectra are consistent with a single 
power-law shape with the same index 1.16 (Sect. [331 Fig.\^. 

3 - After ~ 120 days the strong soft outburst and the modest 
enhancement of the hard non-thermal component both vanish 
(Sect.\M Figs.\^&[M. 

7. Discussion 

The high-B-field radio-quiet X-ray pulsar PSR J1846-0258 ap- 
pears to be a unique object. Ever since the discovery of high- 
energy X-ray and gamma-ray emission from rotation-powered 
(radio) pulsars, they have been known as stable high-energy 
emitters. PSR J 1846-0258 is the first rotation-powered pulsar 
exhibiting a m agnetar-like radiative outburst together wi t h shor t 
< 0.1s bursts dKumar & Safi-Harbl l2008t iGavriil et all l2008l) . 
This schizophrenic behaviour triggered the discussion whether 
high-B-field rotation-powered pulsars represent a transition be- 
tween normal, lower-field pulsars and magnetars. In this discus- 
sion we will first address its manifestation as a genuine rotation- 
powered pulsar and whether its high-energy characteristics are 
commensurate with those of other young rotation-powered pul- 
sars and theoretical scenarios discussed in literature for the 
production of non-thermal emission in pulsar magnetospheres. 
Then, we will discuss how the characteristics which we deter- 
mined for its magnetar-like outburst point to a likely production 
site and scenarios in the dipole geometry discussed in models for 
rotation-powered pulsars. 

7.1. PSR J1 846-0258 as a rotation-powered pulsar 

Before its brightenin g and since its discovery in X-rays by 
iGotthelfetall ( [2000l) . PSR J 1846-0258 manifested itself for 
many years as a stable pulsar emitting X-rays with ener- 
gies above 0.5 keV. The total spectrum (Kes 75 plus the 
PSRJ1846-0258/PWN system) measured with INTEGRAL 
above 20 keV has a photon index F = 1.80 + 0.06 and an X-ray 
luminosity (20-300 keV) of Lf;J^ = 6.3(3) x 10^^ ■ t; ■ erg/s. 
This luminosity amounts ~ 9.7% of the available spin-down lu- 
minosity assuming isotrope emission (rj = 4:7r) and a distance 
of 10 kpc (dio = 1) and is of the same magnitude, but slightly 
higher than what has been measured for other young pulsars such 
as the Crab (3.0%; adopting a distance of 2 kpc and the high- 
ener gy model, a broken power-law, given in lJourdain & Roquesl 
I2OO8 ) and PSR B1509-58 (5.7%; adopting a distance of 5.2 kpc) 
in the same energy band. Also for energies below 10 keV the 
total luminosity of PSR J 184 6-0258 i s similar to that of other 
young pulsars as discussed bv lSu et al] (I2OO8). Obviously, adop- 
tion of a sig nificantly diffe r ent dis tance to PSR J 1846-0258 (e.g. 
19 kpc by IMcBride et all |2008|) will lead to different conclu- 
sions. 

The time-averaged non-thermal pulsed spectrum has a pho- 
ton index r = 1.20 ± 0.01 (3-300 keV, this work) and the X- 
ray luminosity (2-100 keV) of the pulsed emission amounts 
^2-100 ^ 1.90(2) X 10^"* ■ T] ■ t/^„ erg/s, which translates in a 
spin-down luminosity efficiency of 0.23% assuming emission 
into one steradian (77 = 1) and a distance of 10 kpc (c/io = 1). 
These values are in the wide range measured for other young 
rotation powered pulsars emitting at hard X-rays e.g. Crab, PSR 
1509-58 and PSR B0540-69. For example, PSR J 1846-0258 can 
be compared with the young radio pulsar PSR B1509-58, which 
also harbours a strong polar surface magnetic field with strength 
of ~ 1.5 X 10'^ G, in this case below the quantum critical field 
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Fig. 11. The high-energy pulsed specti-a of PSR J1846-0258 
from pre-outburst observations with RXTE and INTEGRAL 
(green data points; green curved line) and from observations per- 
formed during the early phase of the outburst (Al) from RXTE 
PCA data (purple data points). For comparison purposes the 
pulsed spectrum of the young rotation-powered high-B field pul- 
sar PSR B1509-58 is superposed (red data points from RXTE 
PCA and HEXTE, INTEGRAL ISGRl and CGRO BATSE and 
COMPTEL; red curved line, fit to RXTE PCA and HEXTE data) 
along with the pulsed spectra (RXTE and INTEGRAL) of AXPs 
4U 0142+61 (orange data points) and IRXS J1708-4009 (aqua 
data points). Error bars are Icr. All spectra are corrected for in- 
terstellar absorption (unabsorbed). 



strength. Furthermore, PSR B 1509-58 has a similar broad single- 
pulse pro file at soft and hard X-rays, and is detected up to at least 
10 MeV (iKuiper et al.L Il999t ICusumano et al.L I200I . In the 3- 
300 ke V energy band, the shapes of the pulsed non-thermal spec- 
tra of both pulsars are similarly curved (see Fig.fTTT). 

The main difference is that the pulsed flux of 
PSR J 1846-0258 is about a factor of ten lower than the 
pulsed flux from PSR B1509-58 (see Fig. [IB- Taking the 
different source distance s into account, 5.2 + 1.4 kpc for PSR 
B1509-58 and ~ 10 kpc (ISu et all lloosh for PSR J 1846-0258, 
the luminosities differ only by a factor of ~ 2.7. However, 
the efficiencies to convert spin-down energy into X-rays 
are (very) comparable, ei846/ei5{)9 ~ 0.8 (the spin-down 
luminosity of PSR B 1509-58 amounts 1.77 x lO" erg/s, 
and that of PSR J1846-0258 8.19 x 10^^ erg/s). Therefore, 
PSR J 1846-0258 does not appear to be intrinsicly different 
from other young normal rotation-powered (radio) pulsars. 

The scenario responsible for the production of non-thermal 
X-ray and gamma-ray emission from rotation-powered (ra- 
dio) pulsars has been under debate for decades. There is 
still no consensus on the origin of the high-energy electrons 
and/or positrons that are responsible for the high-energy pho- 
tons, and how and where they are accelerated in the open 
zone along the magnetic-field-aligned electric field. For many 
years the discussion was between polar-cap scenarios in which 
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the acceleration was proposed to take place near the neu- 
tron star at the magnetic poles dArons & Scharle mannl 1 19791: 
[Oaughertv & Harding, 1982), and outer gap models placing the 
site of acceleration and emission between the null-charge sur- 
face on which the magnetic field becomes perpendicular to the 
rotation axis and the light cylin der on which the p las ma coro- 
tates with the speed of light ( Cheng et all Il986allbl : iRomaniL 
Ll996; Hirotani & Shibata, 2001). Both approaches encountered 
problems, e.g. the low-latitude polar-cap emission was diffi- 
cult to reconcile with double-peaked or wide pulse profiles and 
the outer-gap scenario faced difficulties in reproducing e.g. the 
shape of the Vela pulse profile. 

For the p olar-cap appro ach the solution was found with 
the slot gap dAronsl Il983h . a narrow bundle of field lines 
bordering the closed-field region in which acceleration and 
pair cascades can occur at altitudes up to several stellar radii 
dMuslimov &'Hardingil2003ll2004lDvks et al.il2()04 ). Slot gaps 
are predicted to exist only in the yo ung and/or fast pulsars. In the 
latest version [Harding et alj d2008l) presented results from a 3D 
model in which they include emission produced by primary elec- 
trons accelerated to high-altitudes in the unscreened electric field 
of the slot gap, as well as by higher-generation electron-positron 
pairs cascading from lower altitudes along field lines interior to 
the slot gap. Curvature, synchrotron and inverse Compton ra- 
diation of both primary electr ons and p airs contribute to the 
total broad-band emission. Hard ing et al .l (2008) show that the 
slot-gap model can reasonably well reproduce the Crab-pulsar 
profiles and (phase-resolved) spectra. The optical to hard X-ray 
component for this young pulsar is produced by synchrotron 
radiation of the unaccelerated pairs flowing along the interior 
field lines of the slot gap to high altitudes. Applying a slot-gap 
model also for the young PSR J 1846-0258 synchrotron radi- 
ation of pairs is then responsible for the stable time-averaged 
pulsed hard-X-ray spectrum (pre- and post-outburst) shown in 
Fig.m The broad single pulse profile of PSR J 1846-0258 can be 
modelled in this scenario for diff'erent viewing directions (angle 
between observer and rotation axis, and angle between rotation- 
and magnetic field ax is) than those modelled for the Crab pulsar 
dHardinget al.Ll2008l) . 

For the traditional outer-gap scenario, in which the non- 
thermal emissio n is produced above the null-charge surface, 
IZhang & ChenS (12000) showed in a 3D model that they could 
reasonably well reproduce for the two young pulsars PSR 
B1509-58 and PSR B0540-69 both their broad pulse pro- 
files (similar to PSR J1846-0258 ) was wefl as their soft 
to hard X-ray spectra. However, in the traditional outer-gap 
model the existence of th e outer gap had been hypothesized. 
iHirotani & Shibatal ([1999) demonstrated in a one-dimensional 
analysis of a pair-production cascade along the magnetic field 
line, that the outer gap does exist as solution of the Maxwell 
and Boltzmann equations. This model was further elabo rated in 
a 2D scenario (Takata et al., 2006; Hirotani, 2006; Taka ta et al.L 
l2007al) and finally the photon propagatio n and pair production 
is solved in the full 3D magnetosphere dTakata et all l2007bl: 
iHirotanil |2008|) . For these self-consistent solutions the gap 
extends past the null-charge surface toward the neutron star. 
[Takata et al. (2007b) and Hirotani (2008) showed that the mod- 
ified outer-gap model can also reproduce the puls e profile and 
phase-r esolved spectra of the Crab. In addition, Takata et"aLl 
d2007bl) reproduced the spectrum and broad structured pulse pro- 
file for the young PSR B0540-69. Like for the slot-gap scenario, 
in the framework of the outer-gap accelerator synchrotron emis- 
sion is responsible for the emission in the X-ray band. This syn- 
chrotron component is dominated by the emissivity of accel- 



erated primary pairs, with a substa n tial con tribution from sec- 
ondary pairs. Notably, i Takata et al.l (l2007b l) show that hard X- 
ray synchrotron emission from secondary pairs below the null- 
charge surface, thus close to the neutron-star surface, contributes 
up to energies above 1 MeV. At MeV energies this latter syn- 
chrotron component even dominates over the synchrotron con- 
tribution from primary pairs. 

In fact, the geometry of the modified outer-gap model and 
that of the slot-gap model are now very similar, but the radiation 
mechanisms are still very different. Both models appear to be 
reasonably capable of reproducing spectra and pulse profiles of 
young pulsars. The X-ray spectral and temporal characteristics 
measured for PSR J 1846-0258 before and after the magnetar- 
like outburst are very similar to those seen for other young pul- 
sars, e.g. PSR 1509-58 and PSR B0540-69, and can also be ex- 
plained with either model. 

7.2. PSRJ1846-0258 exhibiting magnetar-iike beliaviour 

We were able to study the magnetar-like outburst in detail ren- 
dering important directions on what most likely caused the event, 
where the initial trigger happened, and what are likely scenarios 
in the magnetosphere for the production of the outburst emission 
above 3 keV. 

7.2.1. Magnetar-like outburst below 10 keV 

Firstly, the evolution of the spectral shape of the pulsed emis- 
sion over the outburst time interval (see Fig. |8]l clearly shows 
that the major ghtch of size Av/v (2.0 - 4.4) x 10"^ instantly 
triggered the ejection of a soft component with an initial flux 
increase by a factor ~ 5 for energies 2-10 keV. The pulsed spec- 
trum between 3 and 30 keV can be described for the first 125 
days (time interval Al, A2 and B in Fig.|2]and Table[3]) with the 
sum of 1) a new soft component which appeared to soften over 
the burst duration (the photon index evolved from 2.96 + 0.06 
via 3.27 + 0.27 to 4.74 + 0.71, indicative for cooling), and 2) a 
non-thermal component with the same photon index of ~ 1 .2 as 
the pre-outburst pulsed spectrum, but initially enhanced in flux 
by ~ 35%. The new soft component vanished at post-outburst 
epochs and the non-thermal component returned back to its pre- 
outburst flux levels. We do not have spectral information on the 
pul sed emiss ion below 3 keV, but Kumar_& Safi-Harb (200j) 
and iNg et alj (2008) show evidence for a black-body (BB) com- 
ponent in the Chandra total point-source spectrum (pulsar plus 
DC) seven days after the start of the outburst (during interval Al 
in Fig. I2I1. The estimated temperatures were in the range kT ~ 
0.7-1.1 keV with BB radii in the range ~ 0.2-0.6 km, in agree- 
ment with the conventional polar-cap radius. These authors did 
not detect a BB component in 2000, when PSR J 1846-0258 was 
in its low persistent state. We do not consider this non-detection 
as being due to low statistics, because the high-energy X-ray 
and gamma-ray spectra of young pulsars are known to be dom- 
inated by non-thermal emission and do not exhibit a BB en- 
hancement. The latter becomes apparent for middle-aged pulsars 
like Vela and older pulsars. An exception is PSR Jl 119-6127, 
which is also a young high-magnetic-field pulsar with spin prop- 
erties very similar to those of PSR J1846--C)258, but with differ- 
ent X-ray properties. ISafi- Harb & Kuma^ d2008l) show that the 
spectrum of PSR Jl 119-6127 has a thermal component which 
can be fitted by either a BB model with a temperature kT ~ 
0.21 keV, or a neutron star atmospheric model with temperature 
kT ~ 0.14 keV. Both temperatures are significantly lower than 
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that measured during outburst for PSR J1846-0258. Therefore, 
we conclude that the BB component of PSR J1846-0258 dur- 
ing outburst is like the soft component measured with RXTE 
initiated by the major timing glitch and cools down in ~ 120 
days. Then, the measured soft transient component of the pulsed 
emission between 3 and 10 keV in excess to the underlying 
hard non-thermal component can be due to resonant cyclotron 
scattering (RCS) thereof, as has been propose d for e.g. magne- 
tars for their persistent and transient emi ssions dThom pson et al., 
l2002VBelob orodov & Thompsonll2007l) . For example. ,Rea et ah 
(2008) applied this model to the total X-ray spectra of 10 magne- 
tars, including canonical and transient AXPs and SGRs. Fig.fTTI 
shows also the pulsed high-energy spectra of AXPs 4U 0142+61 
and IRXS J 1 70849-4009 10, again m easured with RXTE and 
INTEGRAL dden Hartog et all l2008a''P), to be compared with 
the pulsed spectrum of PSR J1846-0258 directly after the ma- 
jor glitch. The similarity in shape is evident. The pulsed spec- 
tra between 3 and 30 keV of 4U 0142-H61 and IRXS J170849- 
400910 can also be represented by the sum of two, a soft and 
a hard, power-law components. The soft component can plausi- 
bly be explained with a resonant cyclotron scattering origin. The 
hard non-thermal power-law spectra of 4U 0142+61 and IRXS 
J 170849-4009 10 have been measured with INTEGRAL up to 
~ 200 keV and could have a similar origin in the outer mag- 
netosphere as the hard non-thermal X-ray spectra measured for 
young r adio pulsars, as has been discussed in den Hartog et al. 
( l2008bl) . Furthermore, and most interestingly, we found that dur- 
ing the outburst the pulse profile in this energy range below ~ 9 
keV, while completely dominated by the soft flux enhancement, 
is identical to the pulse profile before and after the outburst when 
the emission was purely non-thermal (Fig. |5]l. This means that 
the soft enhancement is produced in the same viewing direc- 
tion through the pulsar magnetosphere. Unfortunately, there is 
no measurement of the profile shape during outburst below 3 
keV, where the soft component dominates. 



7.2.2. Flux enhancement above 10 keV 



In Section l54l we showed that above 10 keV, the evolution over 
the outburst period of the pulsed spectrum differs from that be- 
low 10 keV. Namely, the pulsed flux (10-30 keV) increased af- 
ter the glitch by only ~ 35%, returning to the pre-outburst level 
in ~ 120 days together with the quenching of the new soft 
component. The pulsed spectrum (10-30 keV) remained hard 
during the outburst, consistent with the pre- and post-outburst 
index value of ~ 1.2. Unfortunately, INTEGRAL did not ob- 
serve PSR J1846-0258 when the glitch triggered the outburst. 
Nevertheless, INTEGRAL observations during the tail of the 
outburst covered also PCA observation-time window B (Table 
[3] and see overlap in Fig. |7]l, but the statistics were insufficient 
for extracting the pulsed signal. They revealed, however, an in- 
crease of the total flux (20-100 keV) with (52 + 23)% and no 
change in spectral index. This is a indication that the hard X-ray 
spectrum from PSR J 1846-0258 changed only in normalization 
and not in shape over the broad 10-100 keV energy window, 
suggesting that the production mechanism of the non-thermal 
emission in the pulsar magnetosphere did not change during the 
outburst. Furthermore, the shape of the pulse profile above 10 
keV, as measured with the PCA, also did not change, suggest- 
ing that the geometry of the production site of the non-thermal 
emission in the pulsar magnetosphere did not alter either 



7.2.3. Production scenarios during magnetar-like outburst 

Summarizing, our results strongly suggest that during the out- 
burst, there are two pulsed components: 1) a transient soft com- 
ponent which dominates the pulsed emission below 10 keV, and 
exhibits the same pulse profile as the pre-outburst non-thermal 
emission; 2) a non-thermal hard component with the same char- 
acteristics (spectral shape and pulse profile) as the pre- and post- 
outburst pulsed emission, but with an ~ 35-% enhanced flux av- 
eraged over 32 days after the major glitch. Apparently, the ge- 
ometry of the magnetic field in the magnetosphere did not vary, 
i.e. the assumed pulsar dipole field did not change during the 
magnetar-like outburst. 

During the outburst, the pulsed spectrum of 
PSR J 1846-0258 looks very much like the spectra of the 
persistent emissions from AXPs, as we showed in Fig. [TT| 
for AXP 4U 0142+61 and IRXS J 170849-4009 10. For AXPs 
strong variability and outbursts were detected for energies below 
~ 10 keV, while the hard non-thermal power-law emissions 
measured with INTEGRAL up to ~ 200 keV remained stable 
within ~ 20% on time scales of 0.5-1 year (den Hartog et al.L 
2008a b). Attempts to see correlated flux variability between 
the soft emission below 10 keV and the hard emission above 
10 keV on shorter time scales were not succesful due to the 
poor statistics in all observations above 10 keV with RXTE 
HEXTE and INTEGRAL ISGRI. Therefore, similar coiTelated 
flux variability for AXPs can not be excluded. The abrupt flux 
enhancements of AXPs (and SGRs) below ~ 10 keV have been 
interpreted as a sudden release of energy, either from below the 
crust or above. Apparently, we witness here for the first time 
a similar event in the case of a high-B-field rotation-powered 
pulsar, and the following scenarios emerge, which can explain 
the soft-hard-correlated flux variability. Namely, sufficient stress 
has been built up in the crust to cause cracking, resulting in a 
major glitch, releasing at the same time a significant amount of 
energy, generating the BB -component measured with Chandra 
in ffie total emission during the outburst. Resonant cyclotron 
upscattering thereof subsequently generates the decaying / 
cooling soft pulsed component measured below 10 keV at the 
phase of the pulse of the non-thermal component. This suggests 
that the event causing the abrupt rise, occured near the polar cap 
(no change in pulse profile). In the two models discussed above 
for rotation-powered pulsars, slot-gap and out-gap models, the 
generation of the pulsed non-thermal emission starts above 
the polar cap near the neutron star surface. The modest flux 
increase of the non-thermal emission can then be explained 
in both model scenarios: 1) In the case of the slot-gap model, 
an abrupt increase in the number of electron-positron pairs 
above the polar cap, will temporaly increase the number of 
unaccelerated pairs flowing along the interior field lines of the 
slot gap to high altitudes producing the synchrotron radiation 
responsible for the non-thermal hard X-rays. This provides a 
natural explanation for the stable pulse profile, because the 
geometry did not change. 2) For the extended-outer-gap model 
we are similarly dealing with enhanced synchrotron emission, 
namely from the abrupt enhancement of secondary pairs created 
above the polar cap below the null-charge surface, which are 
accelerated in the extended gap and contribute up to energies 
above 1 MeV. 

The above scenarios would imply ffiat PSR J1846-0258 did 
not transform from a rotation-powered pulsar into a magne- 
tar, rather that rotation-powered pulsars with a strong B field 
can occasionally experience major glitches, releasing signifi- 
cant amounts of energy from below the crusts, following subse- 
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quently one of the above sketched scenarios. With the increase Zhang, L. & Cheng, K.S. 2000, A&A, 363, 575 
of the B-field strengths going from rotation-powered pulsars to 
AXPs and SGRs, the frequency of such events increases. 
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